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Abstract 

At finite temperature we study the quantum tunneling of magnetization for 
a small ferromagnetic particle with the biaxial symmetry placed in a mag- 
netic field at an arbitrary angle. We present numerical WKB exponent below 
the crossover temperature in which the quantum tunneling is affected by the 
thermal activation, and the approximate form of the WKB exponent around 
the crossover region. The effect of quantum fluctuations on the thermal acti- 
vation rate beyond the crossover regime are discussed. 
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I. INTRODUCTION 



In the last decade, the problem of quantum tunneling of magnetization(QTM) in a single 
domain magnetic particle has attracted a great deal of theoretical Jl[ and experimental |2| 
interest. There are two reasons for this. One is the rapid development of the new field of 
nanomagnetism, in which a single domain magnetic particle can be prepared. The other is 
that such a system meets the main criteria for observing macroscopic quantum tunneling 
(MQT). There are three conditions for seeing MQT which can be stated qualitatively as 
follows. First, the height and the width of the barrier should not be large. Second, the 
effective moment of inertia associated with the magnetization switching should not be too 
large. These two conditions are required for a tunneling time accessible to current exper- 
imental situations. Third, the crossover temperature at which the transition between the 
thermal activation and the quantum tunneling occurs should be of the order of milliKelvin 
range or above. Based on the criteria, it has been believed that a magnetic field is a good 
external parameter to make QTM observable. Also, tunneling dynamics poses many inter- 
esting problems, particularly concerning the role played by the thermal fluctuations at finite 
temperature. This is the main subject of this paper, which deals with thermal corrections 
to quantum tunneling rate of a single domain ferromagnetic paricle in a magnetic field. 

Consider a single domain ferromagnetic particle uniformly magnetized along an easy axis 
determined by magnetocrystalline anisotropy. Since the anisotropy generates two or more 
energetically equivalent orientations of the magnetization, at extremely low temperature 
the direction of the magnetization M might change from one easy axis to the other easy 
axis by the quantum tunneling process. However, according to the criteria, the rate of 
change of M is too small to be observed without controlling the height and the width of 
the barrier as well as the effective moment of inertia intrinsically produced by the magnetic 
anisotropy energy. Applying an external magnetic field H in a proper direction, one of two 
equivalent orientations becomes metastable and the magnitude of the physical quantities can 
be controlled in an appropriate way. In this situation the change of M with time is expected 
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to be able to be observed and compared with theroretical results. As temerature increases 
from zero, thermal effects add to the quantum tunneling process. H At sufficiently high 
temperature, the direction of the magnetization is changed by the pure thermal activation 
whose rate is proportional to exp(— U/ksT) where U is an energy barrier related to the 
magnetic anisotropy energy and the external magnetic field. In the intermediate temperature 
region, it is expected that there exists a crossover temperature T c . Around T c both thermal 
fluctuations and quantum tunneling coexist, namely, either the thermally activated quantum 
tunneling below T c or the quantum mechanically assisted thermal activation process above 
T c can occur. Well below the crossover temperature, i.e., at a temperature low enough to 
neglect the thermal activation process, the rate of change of M is soley due to quantum 
tunneling. [|,[|] 

A number of theoretical studies about QTM for the single domain ferromagnetic particle 
have been performed for several magnetocrystalline anisotropy. The QTM problem with 
a magnetic field at an arbitrary angle was firstly studied by Zaslavskii || who calculated the 
tunneling rate for the uniaxial symmetry by mapping the spin system onto a one-dimensional 
particle system. For the same symmetry, Miguel and Chudnovsky [[| have calculated the 
tunneling rate within the imaginary time path integral method. J7| They have also discussed 
the tunneling rate at finite temperature and suggested experimental procedures. Kim and 
Hwang || have performed the calculation based on instanton approach for biaxial and 
tetragonal symmetry. Their work presented the tunneling and oscillation rate between angles 
tt/2 < 6 H < 7r at zero temperature, where 9jj is the angle between the initial easy axis and 
the external magnetic field. This paper extends the tunneling rate for the biaxial symmetry 
to a finite temperature. Within the instanton approach, we present the numerical results 
for the WKB exponent below the crossover temperature and their approximate formulas 
around the crossover temperature. Also, we discuss the effect of quantum fluctuations on 
the thermal activation rate beyond the crossover regime. 

This paper is structured in the following way. In Sec. [Q|, we introduce general formulation 
for the tunneling rate based on the spin coherent state path integral method. We discuss the 
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angular dependence of the critical field and the critical angle by using approximate formula 
of the total energy in the small e(= 1 — H/H c ) limit, where H c is defined as a critical 
field at which the barrier disappears. In Sec. [[II], we give a brief survey of the quantum 
tunneling of magnetization in the presence of the magnetic field with a general direction at 
zero temperature. In Sec. [TV], we consider thermal effects on the quantum tunneling rates 
below T c and qunatum corrections to thermal activation rates above T c and present their 
numerical results in the entire temperature regime. In Sec. |V| we summarize the analytic 
results discussed and provide actual estimates of such important quantities as the tunneling 
rate, crossover temperature, and so on, for several real materials. 



II. GENERAL FORMULATION 

Consider the spin coherent state path integral representation of the partition function 
given by 

Z((3h) = f D[M(r)]exp(-S E /h), (1) 

where ft = l/k B T, the path sum is over all periodic paths M(r) = M(r + fth), and Se the 
Euclidean action [[!]] which includes the Euclidean version of the magnetic Lagrangian L E as 

S e [M(t)] = fdrL E [M(r)]. (2) 

Due to strong exchange interaction in a single domain ferromagnetic particle, the magni- 
tude of the magnetization is a constant M . For that reason the dynamical variable is the 
direction of the magnetization, equivalently, 8(t) and 0(r) in the spherical coordinates of 
M. Therefore, up to the normalization the functional measure in Eq. (|l]) is equivalent to 
DQ as 

N 9 / + 1 

Dn = Jim J] (— ) sin 6 k d6 k d<\> k , (3) 

k=l 

where e = maxfYfc+i — r k ) and J = MqV/Ti^. Here 7 = gfiB/fr, with g being the g-factor, 
fiB the Bohr magneton, M the magnitude of the magnetization and V the volume of the 



system. Since the tunneling rate T of a metastable state is proportional to the imaginary 
part of the free energy of the system, ImF, and ImZ oc ImF, from Eq. ([[]) the tunneling 
rate in semiclassical limit, with an exponential accuracy, is 

r cx exp(-^ in (T)/^), (4) 

where S^ m (T) is obtained along the trajectory with period f3Ti that minimizes the Euclidean 
action 

S E = V f n {i^[l- CO sB{T)\^P-+E[M.{r)\}dT. (5) 
jo 7 cLt 

The first term in Eq. (|^) is the topological Wess-Zumino term |||9|], and the second term is 

the energy density, which is composed of the magnetic anisotropy energy E a and the energy 

given by an external magnetic field H, given by 

E[M(r)] = E a - M ■ H. (6) 

We have selected the biaxial symmetry with magnetic anisotropy energy E a = Ki(a\ + 
a 2 ,) + K 2 a 2 . Here a^s are the directional cosines and K\ > 0, K 2 > are the parallel and 
transverse anisotropy constants, respectively, whose relative magnitude can be of any value. 
Applying the magnetic field in the xz plane [|,j5| and expressing «i, a 2 and «3 in spherical 
coordinates of M, the energy (^) is expressed as 

E[6, 0] = K x sin 2 9 + K 2 sin 2 sin 2 9 

-H X M sin9 cos - H Z M cos 9 + E , (7) 

where E is a constant to make E(9, 0) zero at the initial orientation. As will be seen 
later, while there isn't an exact analog of kinetic energy in the action (||), there is an 
effective moment of inertia in the dynamics of a single domain particle that is inversely 
proportional to a linear combination of K 2 and H x . For this reason we need either the 
transverse anisotropy constant K 2 or the magnetic field H x transverse to the initial easy 
axis for quantum tunneling. Otherwise, QTM is not observable due to the effective moment 
of inertia whose magnitude becomes infinity. 
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Introducing the dimensionless constants, 



k 2 = K 2 /K u h x = H x /H , h z = H z /H , H = 2K 1 /M ( 



(8) 



we obtain the energy ([?D written as 




— h x sin 9 cos 4> — h z cos 9 + Eq, 



(9) 



where E(9, 0) = E(9, <f>)/2K\. As noted in Eq. @, there are two equivalent easy directions 
in the xy plane, which are = and 7r in the azimuthal angle if h x = 0. We treat the 
problem of the easy directions given by = in the xy plane. Even though we start the 
argument from = ir, the whole discussion is the same as the one from = by replacing 
by 7T + 0. So, the classical path becomes c i(t) = 0o( r ) + nn where n = or 1 where 0o(t) 
is the classical trajectory started from = 0. || However, if the external field is not along 
the easy axis(z axis), i.e., h x ^ 0, only = in Eq. (||) is the easy direction in the xy plane 
and = 7r not. Keeping this point in mind, the total energy on the easy plane becomes 



where we used h x = h sin 6 H and h z = h cos 9h- Let us define 9 to be the angle of metastable 
state generated by the anisotropy energy and the external magnetic field, and 6 C the an- 
gle at which the barrier vanishes by the external critical magnetic field H c . Then, 6q is 
determined by [dE(9,0)/d9]e=e = 0, and 9 C and h c by both [dE(9, 0)/d9]e = o C! h=h c = 
and [d 2 E{9 1 Q)/d9 2 }e=e c ,h=h c = 0. After simple calculations, the dimensionless switching or 
critical field h c (= H c / H ) and the critical angle 9 C are expressed as [UH 




(10) 



/i c = (sin 2 / 3 ^ + |cos^| 2/3 )" 3/2 
1 . r 2\cot9 H \ 1/3 , 




(11) 



(12) 



where it should be noted that Eq. flTT| ) is the well-known Stoner-Wohlfarth expression. 
|T0| Simple analysis for Eq. ([H]) reveals that the small value of e is favorable for the 
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magnetization switching because the height and the width of the barrier are proportional to 
the power of e. The situation can be achieved by applying the external magnetic field close 
to the critical field. Hence the problem will be considered in the small limit of e, i.e., e<l. 

Expanding [dE(0, O)/d0]g=g o = about C , and 

using the relations [dE(0,0)/dO]g = g c ^=h c = and [d 2 E(0, 0) / d0 2 ] e=ecjh=hc = 0, we obtain 
the equation for 77 (= 9 C — 9 ) as [|] 

sin(2# c )(e - -rj 2 ) - 77cos(2# c )(2e - r, 2 ) = 0. (13) 

Simple calculations show that 77 is of the order of y/e. Thus the order of magnitude of the 
second term in Eq. ([13]) is smaller than that of the first term by \fe and the value of 77 is 
determined by the first term which leads to 77 ~ y^2e/3. However, when 9h is very close 
to 7r/2 or 7r, sin(2# c ) becomes close to zero, and the second term is much larger than the 
first term. Therefore, the value of 77 is obtained from the second term when 9h — tt/2 or tt, 
which leads to 77 ~ \/2e for Oh — vr/2 and 77 ~ for Oh — vr from the detailed analysis of the 
equations which 0q and C satisfy. Since the first term in (|T3"D is dominant in the range of value 
C which satisfies tan(2# c ) > 0(y/e), r] ~ y / 2e/3 is valid for tt/2 + 0(y/e) < H < vr - 0(y/e) 
from Eq. (|12|). Fig. |1] shows that this is checked by performing the numerical calculation 
for the equations which 0q and C satisfy. 

Introducing a small variable 8(— — 0q), we obtain an approximate form of E(0, 0) as 

E(6,<p) = ysin 2 0sin 2 (e o + 5) + /ixSin(^o + 5)(l-cos0) + J E 1 (5), (14) 

where E(0, 0) is represented as E(S, 0), and Ei(S) is a function of only S given by 

E x {5) = ^sm(20 c )(35 2 7]-5 3 ) 

+\ C o^20 c )[5 2 {e-\r l 2 ) + 5S-^}- (15) 

Although cos(26 l c )-term in Eq. (|I^) looks smaller by a factor of 77 which is of the order of 
v^e, it can not be neglected near Oh = tt/2 and tc because sin(2# c ) is almost zero for these 
regions of Oh- In order to evaluate the order of magnitude of each term in action (|5|) and 
simplify the calculations for e <C 1, it is convenient to use new scaled variables 
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f = e a/2 u T, 5 = 5/y/e, lo = 2 1 K 1 /M . (16) 



Then, the Euclidean action ([5]) becomes 



S E [8(f),<i>(f)] = hJe^ ( df{iet[l - cos(fl + v^)]-? 

Jo ar 

+— sin 2 sin 2 (^o + y/e5) + h x sm(8 + \^5)(1 — cose 



+ \ sin(2^)e§ (3?^= - f) + i cos(2£ c ) e 2 [5 2 (l - ^) + 5 3 ^= - ^]}. (17) 

where (3 = e a ^ 2 uo(3Ti and the parameters a is fixed by the analysis of the order of magnitude 
of each term in the action fllTD subject to the situation considered. 



III. QTM AT ZERO TEMPERATURE 

Consider the Euclidean action ( |TTD at T = || by performing a simple dimensional 
analysis. At T = (J3 — > oo) the instanton starts and ends at the metastable state 5 C \ = 0, 
0d = for tt/2 + 0{y/e) < 9 H < vr - 0{y/e), in which E(6 d = 0, cl = 0) = 0. Thus, from 
the energy conservation the total energy E(5 C \, <p c \) in Eq. ([14]) becomes zero at any f for 
T = 0. Since sin(2# c )-term (~ e 3 / 2 ) is larger than the last term (~ e 2 ) by y/e in Eq. (|17|), 
the first, second and third term should be of the order of e 3 / 2 . Since h x is finite in the range 
of angle and #0 is not close to and 7r/2, the small value of <p in the third term of Eq. 

contributes to the path integral such that <p should be of the order of e 3//4 . This fact 
is valid for any value of because h x is not zero in this range of angle and plays the role 
of the effective transverse anisotropy component which is crucial for QTM. Therefore, from 
E(6 c i, c i) = 0, the magnitude 4> c \ is approximately given by 



9 = 2£ 1 (* cl ) 

rl k 2 sin 2 (9 + S c i) + h x sin(0 o + 5 C 



(18) 



From Eq. ( |15D we get Ei(5) ~ sm(29 c )(3S 2 rj — 5 3 ) because cos(2# c )-term is much smaller 



than sin (29 c ) -term. Since 77 ~ y2e/3, and the critical angle 9 C and the angle 9q of metastable 
state have values between and tt/2 in this range of angle, we get <p c \ ~ — i^/M cot # c e 3//4 , 



where the effective moment of inertia M is 11 



Thus, the approximate formula of the classical Euclidean action is found to be S E ~ 
ih J sin Q c \fe5 c \(j) c \ ~ % Je 5 ' 4 \/M sin C cos C from the dynamical part of the first term in Eq. 
(|l7l), where 5 C \ ~ 0(1) and E(5 C \, c i) = 0. By using the critical field (|Tl"l) and the critical 
angle (|12"1), the classical action is approximately 



Si « ftJe 5 / 4 ^), (20) 

where 

s(M = „ (21) 



1 + g (1 + | COt^| 2 /3) 

In the case of Oh = 7r, we have 6* c = 6* = ?7 = and < 5 C \ < 2y/e from Eq. (|T5|). In 
this situation, the action ( [FTP becomes 

S,[I(f),*(f)] = hJe-i [if{\^f (§) + *J>*?* + ^<F - J)}. (22) 



Since the last term in Eq. fl22[ ) is of the order of e 2 , the order of magnitude of fc 2 sin 2 in 
the second term should be e. If k 2 becomes much smaller than e, in which the magnetic 
particle possesses the uniaxial symmetry, f||| QTM is not observable. This is understood 
from the fact as follows. Consider the action ((5|) as the dynamical system in the Hamiltonian 
formulation which consists of the canonical coordinate and p^ = 1 — cos#. [Q The Poisson 
bracket {p^, H} which determines the dynamics of the spin system with the Lagrangian 
(^) becomes zero for Oh = tt because the Hamiltonian H becomes a function of only p^. 
Thus, does not depend on time. If the value of k 2 is much greater than e, the small 
value of 0, i.e., 0(0) ~ e 1 / 2 mainly contributes to the path integral, in which a is chosen 
to be 1. In this case, a comparison of fc 2 -term with the last term in Eq. (E2f) determines 



0ci ~ —i\jKie/K2, which leads to the least action S E ~ ihJe5 c \ <p c \ ~ %J(?l 2 yjKx/K 2 . The 
approximate forms of S E which we have discussed for n/2 + 0(y / e) < H < n — O(^Je) and 
Oh = 7t agree with the exact results in Ref. || up to the numerical factor iV = 16 x 6 1 ^ 4 /5 
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and 8/3, respectively. |13| From above estimates we note that the dependence of the WKB 
exponent on Oh, e and K 2 /Ki can be deduced from the simple dimensional analysis based on 
the energy conservation E(5 C \, <fi c \) = without solving the equation of motion from SSe = 0. 
However, in order to obtain the dependence of the WKB exponent on K\, K 2 , e, and Oh 
at finite temperature, we cannot make use of the dimensional estimate as above, because 
E(8 C \, c i) 7^ 0. Thus, we need to calculate the bounce solution by numerically solving the 
equation of motion satisfied by the least trajectory 

IV. QTM AT FINITE TEMPERATURE 

Consider the Euclidean action (|17D with a periodic instanton for n/2 + 0{\fe) < Oh < 
it — 0(y/e) and Oh = 7r. 

A. tt/2 + O(Ve) <9 h <tt- 0(y/e) 

Since the sin(26 l c )-term is much larger than cos(20 c )-term in Eq. ([T7|), the action is 
simplified as 

S E [5(f), 0(f)] = ftJe-f / df{-it^ sin(^ + (— ) 

JO CLT 

+ y sin 2 sin 2 (0 O + v^5) + 2h x sin(0 o + y/e5) sin 2 (^) 

+ ]sm(20 c )e-2(3f^-f)}, (23) 
4 ye 

where h x ^ and /C2 has any value, and we performed the integration by part for the first 
term. Also, the total time derivative is neglected because c i(— (3/2) = cl (/3/2) = by 
noting that from the Euler-Lagrange equation c i is proportional to dS c \/df which is zero at 
f = ±/3/2 due to the periodicity of 5 c i(f). As is well known for the crystal symmetry such 
as the form K z a 2 z — K y a 2 , the topological term d^jdr in Eq. ([5]) plays important roles in 
the presence of a clockwise or counterclockwise winding over the barrier along the passage 
from cf) = to = 7r in macroscopic quantum coherence(MQC). Since there is no such 



10 



topological situation in our quantum tunneling problem, the total derivative term does not 
contribute to the classical action. In the action ( |23| ) the last term generates the potential 
barrier in the problem. Thus, the magnitude of each term should be the same order as that 
of the last term. Noting that the order of magnitude of each term is 0(e^ a+1 ^ 2 (p), 0(sin 2 0), 
O(sin 2 (0/2)) and e 3 / 2 , respectively, the small values of <p contribute to the path integral, 
and its order of magnitude is expected to be e 3//4 , as was discussed at T = 0. Therefore, an 
appropriate choice of a is 1/2, which makes to each term to be of the order of e 3//2 in the 
integrand. 

Performing the Gaussian integration over <p in the partition function ([!]) and the measure 



and introducing the variable r = f^Jsm(28 c )/M, the remaining integral is of the form 

jD[5(f)}eM-Sf/h), (24) 



where (3 = /3ysin(20 c )/M = e 1 ^ A uJo/3h^sm(29 c ) / M and the effective action is given by 

Sf[S(f)] = hJe 5 /^Msm(29 c ) j* df [\{^] 2 + U(6)), (25) 

where U(8) = (V6<f - <f )/4. 

The classical trajectory S c \(f) with period (3 which minimizes the effective action satisfies 

d 2 S d r dU, 



df 2 L d5 

and the corresponding total energy becomes 



(26) 



E to m = \(VG?-?)-\(§) 2 . (27) 

At zero temperature the classical trajectory becomes a regular bounce S c \(f) = 
v^6/ cosh 2 [(^) 1//4 f] with E tot ((3 — > oo) = 0. In order to include thermal corrections to 
quantum tunneling rate, we need to consider the periodic solutions which satisfy Eq. (p6|). 
In fact, there exist two kinds of periodic solution, 5 m (= 2y/E/3) and the periodic motion 
^ci(^) of the particle in the inversed potential with total energy E tot , as is seen in Fig. 2. 
|j~4i For the constant solution S m , we obtain the classical action from Eq. (|25|) 
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S c i = S = hJe 5/4 sjMsin26 c (3U{8 m ) = /3hU m , (28) 

where 

^ = ~ ^ lVe [ iTw ]l (29) 

and the corresponding escape rate 

T oc exp(-S /K) = exp(-U m /k B T), (30) 

which is the Boltzmann formula representing a pure thermal activation. Let us define T c to 
be the temperature at which the periodic solution 5 c \(f) with period (3 approaches the limit 
5 m . Then, slightly below T c , the thermal bounce S c \(f) reduces to small oscillations near the 
bottom of the inversed potential —U(5) shown in Fig. |2|. In this case we take only the first 
Fourier harmonics for the thermal bounce because the next harmonics are smaller near T c , 

m 



5 cl (f) = 5 (T) + 5 1 (T)cos(ujf), (31) 

with Co = 2tt/P. The temperature dependence of and 5± is crucial for both the magnitude 
of T c and the dependence of S c \ on f . If S± does not depend on T, 6 C \ becomes a function 
of f through cos(ci} c f ) at T = T c , which is contrary to the fact that 5 C \ is a constant 5 m 
at T c . Therefore, 8\ should vanish at T c and have a finite value below T c . In order to 
obtain T c and the bounce solution for T < T c , we substitute the harmonic type solution 
(j31|) into the equation of motion for the thermal bounce, such as Eqs ([34]) or (f£8) subject to 



the situation under consideration. We then obtain the temperature dependence of 5q and 
5 1 by taking the constant term and the coefficient of cos (cur) to be zero, and the relation 
u c {= 2tc/(3 c ) = -U"{8 m = 2V6/3) which gives 

k B T c =^f(e,9 H ), (32) 

7T 

where 
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,3 u/4a/4r |cot^r I K 2 _ 





cot 9h 


1/6 


1 + 


cotO 




2/3 



/Mi/Hu^v a i , 2/3 ]^i+^(i+icot6) g p/3). (33) 



Since the thermal bounce ( pT|) satisfies the equation of motion 



the temperature dependence of <5 and 5i is given by 



5o(T) = ^[l + (|) 2 ], - Sl (T) = ^l-(^)*, (35) 

where at T = T c 5(f) = 2^/3 (= <5 m ) by noting that 8 (T C ) = 2^/6/3 and 5i(T c ) = 0. 
Substituting in Eq. (^) for the action, the approximate form of the minimal action can be 
written as 

Cmin TT rp rp 

-ir~b [l - 3( -^ fl (36) 

It is important that for T < T c the action S^ m of the thermal bounce ([H]) is smaller than 
the action Um/ksT of the constant path 5 m . Hence below T c , the functional integral for 
the decay is dominated by the thermal bounce. This fact implies that T c is the crossover 
temperature from the quantum-mechanical to the thermally activated decay Noting that 
the Boltzmann formula (|30|) is derived from the constant path S c \(f) = 5 m and also valid 
above T c , the thermal bounce is degenerate into the constant trajectory for T > T c . 

In order to obtain the WKB exponent in the entire range of temperatures less than T c , 
we employ two equivalent approaches for the numerical calculation of the least action. First, 
we use the equation of motion (|3"4| ) for the numerical bounce solution, in which the action 
has been integrated over with a periodic boundary condition 8 c i(f = —(3/2) = 6 c i(f = 
(3/2) where j3 = 27r(2/3) 1//4 (T c /T). As illustrated in Fig. note that numerical solutions 



become oscillatory formulas which are of the form ([U]) at temperatures slightly less than 
T c . When inserted into the action fl25D , the classical action can be obtained by performing 
the integration numerically. Second, numerically solve the coupled equation of motion of 



5 C \ and 4> c \ derived from the Euler-Lagrange equation for the action ( p3|) by incorporating 
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the energy conservation E(S C \, 4> c \) = 0. The second method is reduced to the first one 
if |0| is small. However, since the first method is not valid unless \<p\ is small enough to 
perform the Gaussian integration over <ft, the second method is a more general approach to 
obtain the WKB exponent. Although the first method is enough for the WKB exponent in 
our cases, the second method is useful for checking the results of the first. Thus, we used 
both approaches for the WKB exponent and found that the results obtained from the first 
methods are identical to the one from the second method. 

In Fig.fl we present the phase diagram of 5 and 0(= 4>/4>o) where 0o — 



ig(e H )sJ2(l + | cot^| 2/3 )e 3/4 in the range of temperature < T < T c . Note that the 
orbit is oval for < T << T c and circular near T c with a center at (2v^6/3, 0). Substituting 
S c i(f) and (p c i into the action ([23|), we obtain the WKB exponent which is shown in Fig. |5|. 
Note that the approximate form fl36|) is valid at the temperature close to T c , and the WKB 
exponent is not sensitive to T for temperatures between and T c . 

We now discuss the effect of quantum fluctuation on the thermal activation at the tem- 
perature above T c . At temperature slightly above T c quantum correction to the thermal 
activation becomes important. Its effect is due to the extension of the regular second-order 
action by terms up to the fourth order, i.e., non-Gaussian terms in the action functional. The 
quantum fluctuation, therefore, smears the transition in the crossover region |T C — T\ <C T c , 
in which the WKB exponent is of the form (|36|). Even at higher temperature where thermal 
activation prevails, the quantum effect is still incorporated into the preexponential factor of 



the thermal activation rate in the form |T7 

UJ 



T = -^C q exp(-(3U m ) (37) 

where uj p = \JU"{5 = 0) is the frequency of oscillation around the metastable minimum and 
C q a quantum mechanical correction. The escape rate is enhanced by quantum fluctuations. 
This correction has been studied in great detail [16,n]. Calculation of C p based on quantum 
fluctuation, the decay rate which includes the quantum effect in this region has the form 

r = £ft;^exp(-/^ m ), (38) 
14 



where u n = 2im/ (3h and Ub(= y —U"(8 m )) is the barrier frequency. This barrier frequency 
characterizes the width of the parabolic top of the barrier. Using the identity 

1 1 00 1 1 2 

t = _L TT U n ,aa) 

2sinh(/5^/2) (5huj ^uj 2 + u 2 J 1 ; 

the decay rate can be expressed as [ ]TT| . [TB1 | 

r=( ^ smh^) 
Z7T sm(prt,cc;f ) /2) 

Since t7"(0) = — 0""(5 m ) in this range of angles, we obtain cj b = u; p = 2u f(e, Oh)- At 
temperatures slightly above the crossover region T — T c <ti T c , we expand sin(f3huJb/2) in 
Eq. (pE0| ) about T c . Noting that Pfi^b = 27r from the crossover temperature (0), we have 
sm(hu b /2k B T) 7r(T — T c )/T + • • • and from Eq. ( |4T)D the escape rate is given by 



T « ^/(e,^)-^— smh[(3hf(e,6 H )uj }exp(-f3U m ). (41) 

7T Z 1 — ±r 



At temperatures well below the pure thermal activation regime, but well beyond the 
crossover region T — T c <C T c , we need to take the escape rate fl40|) as an exponential 
of a sum of logarithms and expand each logarithm in power of f3h. Then, we obtain the 
approximate form of the escape rate given by 

r « ^f(e, 6 H ) exp{-[^% - ^^(|^) 2 ]}. (42) 

7T Kgl 6 KbJ- 

As noted in Eq. (f|^), the quantum effect which is inversely proportional to T 2 , is incorpo- 
rated into the Boltzmann formula {(3U m ) in the exponent which is indicative of quantum- 
mechanically assisted thermal activation process. 

The expression (|4*0|) which includes the quantum fluctuation effects is valid in the entire 
thermal activation regime beyond the crossover region T — T c ^ T c . Its approximate forms 
(fn|) and ([|2|) are valid slightly and well above the crossover temperature, respectively. For a 
pure thermal activation regime in which f3 becomes much small, we obtain the approximate 
relation smh(f3Tiuj p /2) « sm(/3huJb/2). In this case the escape rate is given by 

r = exp(-pU m ), (43) 

7T 
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which is just the Boltzmann formula for the pure thermal activation process from Eq. (|40|). 

Fig. ^| illustrates that we superpose the approximate expressions for the tunneling rates 
on the numerical results and show the region of the validity of the analytic results. Eq. ([42]) 
is applicable in a wide range of temperatures above T c but not appropriate as T approaches 
the crossover region which includes T c . Close to the crossover region, the expression ( f4"T| ) 
is valid in a narrow region. In the crossover region Eq. (|4l| ) is not appropriate because 
of singularity. In this case we need to extend the second-order action by terms up to the 
higher order in order to regularize the divergence, as previously noted, which leads to the 



WKB exponent to be of the form (0). |T9[ When compared with the Boltzmann formula 
([43|) based on the pure thermal activation, quantum fluctuations increase as T becomes 
close to T c from the pure thermal regime. Also, note that while the tunneling rate at zero 
temperature is valid in a wide range of temperature below T c , there exist an appreciable 
amount of thermal fluctuations for T 5) T r . 



B. H = 7T 



When the external magnetic field is opposite to the initial orientation, h x = 6 = 9 C = 
and the Euclidean action becomes Eq. (^). Estimating the order of magnitude of each term 
in the integrand of the action (^), a is chosen to be 1 and (f) ~ 0(e 1 ^ 2 ) for the value of k-i 
much larger than e, as previously discussed at T = 0. Performing the Gaussian integration 
over (J), the effective action for the integral (pl|) is given by 

I ,.!<■'■ 2 /•>• ' ' 



(44) 



where we introduced r = \fk2f and (3 = \/fe/3 = e 1 / 2 uJof3hJ K 2 / 'K 1 . At T = 0, the least 



trajectory for the action (03) is 5 c i(r) = 2/cosh(r) with E tot = U(5 C \) - \{^-) 2 = where 



U(6) = {A5 2 - P)/8. 

As with sbusection [IV A| , the height of barrier U m = KiVe 2 , and the crossover tempera- 
ture are 
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, rj, tlUJo I K 2 1/2 

kB% = — \2K x e ■ (45) 



For T slightly less than T c , the thermal bounce is represented as 



Uf) = ^1 + 4(1), + ^1 - (1), cos(^), (46) 



and its corresponding action 

S c e Um fl 32 Tc-T 2 

Ti k B T [ 15 l Tc ' V K } 

In the entire temperature range less than the crossover temperature, the WKB exponent 
is found by using the numerical solution for the equation given by 

;p-« + T = 0. (48) 

with a periodic boundary condition with a period (3 = (2n/y/2)(T c /T). Its result is similar 
to the one illustrated in Fig. £|. Using this solution, we obtain the WKB exponent whose 
form is the same as Fig. ^. 

For the thermal activation regime beyond the crossover region T—T c <C T c , it is necessary 
to consider the escape rate (fit]) by including quantum fluctuation with the frequency uj p = 



€ 1 I 2 loq^I a 2 / Ki(= Ub/y/2). Then, the rate is approximately given by 

T « ^ s (e)-^-smh[(3hs(e)uj ] exp(-(3U m ), (49) 
for the quantum correction regime slightly beyond the crossover region T - T c C T c , and 

r«- a (e)e^{-[^-— (^)]}, (50) 



for the temperature regime well beyond the crossover region, where s(e) = e 1 ' 2 '■ y /K 2 /2K 1 . 
The shape of lnr(T)/r(0) is similar to the results shown in Fig. 19 at all temperatures, 



V. DISCUSSION AND CONCLUSION 



Table | summarizes the analytic results for the range of angles discussed. Note that the 
e-dependence of the height of barrier, the crossover temperature and the WKB exponent at 
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T = are different for tt/2 + 0(y/e) < Oh < it — 0(y/e) and Oh = tt- The exponents of e 
in U m , T c and i?(T = 0) are equal to 1.5, 0.25 and 1.25 in a wide range of angles between 
7r/2 and tt, and increase up to a value of 2, 0.5 and 1.5 if the applied field goes toward the 
opposite direction of the initial magnetization (Oh = tt), respectively. The properties of the 
former are characteristic of the quadratic-plus-cubic potential whose form is well-known in 
Josephson systems, and those of the latter are characteristic of the quadratic-plus-quartic 
potential. [pO 



The exponent of e in the WKB exponent is easily understood by noting that the WKB 
exponent is propotional to the height of barrier U m and inversely propotional to the width 
of the barrier characterized by u;&, i.e., B oc U/hub. To be specific, B oc e 3 / 2-1 / 4 for tt/2 + 
0{y/e) < H < tt - 0(y/e) and B oc e 2_1/2 for H = tt. Denoting the WKB exponent 
B(T = 0) to be CBU m /hu b , cb is 36/5 for the range of angle tt/2 + 0(y/e) < H < n-O(yfe) 
and 32/3 for Oh = tt, which are generally determined by the form of potential near liability. 

The ratio of T c (0h) to T c (135°) is larger at k 2 = than at finite &2 for Oh < 135° and 
smaller for Oh > 135° (Fig. [7]). The position of the maximum of T c moves from Oh — > 101° 
at k 2 — > (uniaxial symmetry) to Oh — > 135° at k 2 — > oo, and its magnitude T c becomes 
smaller as k 2 increases. This behavior is easily understood from Eq. (^) which gives a 
maximaum value of the crossover temperature T c at 



T r — tt — arctan[ , ^ 2/ ^ f\ (51) 

■ / (ir 2 /K 1 + 0.5) 2 + 2-1.5 



as is illustrated in Fig. || 

It is noted that the WKB exponent B at zero temperature can have the maximum 
value, depending on the ratio of K 2 to K\. For K\ ^> K 2 , there is no maximal point in 
B at zero temperature. j|,[5| Thus, there are some values of the ratio of the anisotropy 
constants K 2 /Ki, at which the WKB exponent changes slightly in a wide range of angles. 
In the crossover region the quantum fluctuation is incorporated into the WKB exponent as 
a parabolic form with a slightly different curvature with respect to temperature in the range 
of angle tt/2 + 0(y/e) < H < vr — 0(y/e) and H = tt, which leads to the smearing of the 



WKB exponent, as is seen in Fig. ||. The quantum correction is also reflected in the thermal 
activation regime with different coefficient of (T c /T) 2 well above the crossover temperature 
for tt/2 + 0(y/e) < 9 H < n - 0(y/e) and 9 H = n. 

To illustrate the above approximate results with concrete numbers we have collected 
in Table [TI] and |TTT| various values for several materials, the example with K\ = K2, 
BaFe 12 _2xCo x Ti x Oi 9 (x = 0.8), fn) Tb, [g2j and Ni [f||. From Table |TTl| it is evident that the 
typical magnitude of e for the magnetic system to tunnel out of the barrier within resonable 
time is of the order of 1CT 3 or less, and that the associated crossover temperature T c is 
typically of the order of 100 mK or less. Noting that the inverse of the WKB exponent B^ 1 



is the magnetic viscosity S studied by magnetic relaxation measurements, ||24j the value of 
S(T = 0) at the angle 9h = 135° are of order of 0.1-1 for the parameter e with the magnitude 
10 -3 - 10~ 4 . Also, the estimated values of 6 H C and T c max in Table fn can be compared with 
the experimental measurements and their consistency will suggest the possibility of QTM 
in these particles. 

In conclusion, we have studied thermal effects on quantum tunneling of magnetization 
placed in a magnetic field at an arbitrary angle. We have presented the analytic forms of 
quantum tunneling rates at several temperature regimes and performed the numerical cal- 
culations in the entire temperature regime. It is found that thermal corrections to quantum 
tunneling rate is small for T <C T c , but quite large for T ~ T c . This is because the thermal 
activation process is strongly influenced by the quantum fluctuation in this regime. Fur- 
thermore, the dependence of T c and magnetic viscosity on the parameters 9jj, e and K2/K1 
is expected to be observed in future experiments. 
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TABLES 

TABLE I. Summary of the results for the quantum tunneling in the range of angle 
vr/2 + 0{y/e) < 9 H < tt - 0{y/e) and 6 H = ir. B(= S c ^/K) is the WKB exponent at (a) zero 
temperature (T = 0), || (b) crossover region(T — T c <C T c ), and (c) thermal activation region with 
quantum fluctuations (T-T c > T c ). Here t = T/T c , m(0 H ) = | cot 6 H \ 1/3 /(l + | cot 8 H \ 2 / 3 ), g{6 H ) 
is given by Eq. © and f(e,9 H ) by Eq. ©. 
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TABLE II. Magnetization Mq, easy-axis anisotropy constant Kx, hard-axis anisotropy constant 
K2, critical magnetic field H c , and characteristic frequency ojq for various materials, where in case 
of Ni the shape-induced contribution to the hard axis anisotropy is obtained for the slab geometry. 



M 
[emu /cm 3 ] 



10 5 [erg/cm 3 ] 



K 2 

10 5 [erg/cm 3 ] 



V 

[10 3 nm 3 ] 



HJ135 ) 



[10 3 Oe] 



# c (180°) 



[10 3 Oe] 



K x = K 2 500 

BaFe 12 _ 2x Co x Ti x 3 10 
Oi 9 (x = 0.8) |H 

Tb [||] ~ 100 

Ni H 508 



100 
12 

~ 10 



100 



10 3 



16 



1 

3.14 

~ 1 
~ 1 



20 
3.53 

10 
1.575 



40 
7.06 

20 
3.15 



24 



TABLE III. Angle for the maximal WKB exponent #j| max and the maximal crossover temper- 

j'max 

ature Orf , magnetic field parameter e, the height of barrier U m at Oh = 135° and 9h = 180°, the 
WKB exponent B(T = 0) at 9 H = (9JJ max and 135°, the crossover temperature T c at 9 H = 9 H C , 
and inverse tunneling rate r _1 (T = 0) at 9 H = 135° for various materials. 
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FIGURES 



FIG. 1. r)(= 9 C - 6 ) as a function of 9 H for e = (a) 0.01 and (b) 0.001. Note that 77 » V2e/3 
is valid in a wide range of angles, i.e., n/2 + 0{yfe) < On < tt — 0(y/e). 

FIG. 2. (a) The shape of the potential U(S) and (b) the shape of the inverted potential —U(5) 
for the range of angle tt/2 + 0(i/e) < 9h < vr — 0(-y/e). In (a) the height of barrier U(S m = 2\/6/3) 
is 2^/6/9 at <5 m = 2>/6/3. 



FIG. 3. The bounce solution obtained by numerically solving Eq. (34) for several temperatures 
t = (a) 0.95, (b) 0.98, (c) 0.99 and (d) 0.999, where t = T/T c . The dashed curves in each 
temperature are drawn from the oscillatory solution (pTl) for comparison. Notice that as T becomes 



close to T c , the approximate formulas from Eq. ( |3l|) are compatible with numerical solutions. 

FIG. 4. The phase diagram 4> vs. 6 at t = (a) 0, (b) 0.82, (c) 0.97, and (d) 0.999. Note that as 
T increases from zero, the orbit is changed from an egg-shape to a circle. 

FIG. 5. The temperature dependence of the magnetic viscosity, S(T)(= BiT)^ 1 ): S(T)/S(0) 
versus T/T c . Here the dotted line indicates the analytic formula ( |36| ) in the crossover temperature 
region, \T - T c \ « T c . 

FIG. 6. Plots of lnr/r(0) as a function of T/T c , where T(0) is the tunneling rate at T = 0. (a) 
the numerical calculation in the entire temperature, (b) Eq. ( f43| ) for the pure thermal activation 
regime, (c) Eq. ( f42"D for the quantum corrections, (d) Eq. ( |4l| ) for the quantum corrections close 
to the crossover region, and (e) the rate in the crossover region whose exponent is of the form 



FIG. 7. The 6h dependence of the scaled crossover temperature T c [= T c (8h)/T c (135°)] at (a) 
k 2 = 0, (b) 1.0, (c) 2.0, and (d) 100. 

FIG. 8. The dependence of 9jf on k 2 which is derived from Eq. (51) where k 2 is defined in 
Eq. d). 
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